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CELLULAR AUTOMATON URBAN TRAFFIC MODEL 


Felix TOTIR', Stefan CANTARAGIU”, Tudor POPESCU® 


Rezumat. Articolul abordeaza modelarea traficului urban de automobile. 
Principalele caracteristici ale metodei.expuse sunt utilizarea unui automat celular 
pentru modelarea distributiei autovehiculelor de-a lungul.a doud strazi intre care 
este interpus un semafor. Se analizeaza astfel aparitia si incidenta fenomenelor de 
ambuteiaj. 

Modelarea comportamentului conducatorilor auto se face probabilist. Rezultatele 
prezentate pot fi utilizate atat in vederea asigurarii unei fluente sporite, cat si a 
optimizarii consumului de combustibil, controlul nivelelor de poluare chimica si 
fonica. 


Abstract. The paper presents an approach for urban vehicle traffic modeling. Main 
characteristics of the exposed method are the use of a cellular automaton in order to 
model the vehicle spreading over two streets interconnected with the help of a traffic 
light. The apparition and the incidence of bottlenecks phenomena are thus analyzed. 
Drivers’ behavior is modeled probabilistically. Obtained results are valuable for 
increasing traffic fluency, optimizing fuel consumption and controlling chemical and 
Phonic pollution levels. 
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1. Introduction 


Appropriate modeling of urban traffic is of extreme importance in the modem 
days. Heavily crowded, cities must cope with a growing number of cars, the level 
of traffic continually increasing. Fluency of vehicle traffic is not the only issue: 
levels of phonic and chemical pollution are also to be addressed. Fuel 
consumption is another first-hand issue, especially with the non-renewable, fossil 
resources showing their limits. 


Modeling a city’s vehicle traffic is, quite often, a particular and stand-alone 
problem. This 1s because, for the vast majority of cases, no networks of roads are 
the same for any two cities. To further complicate the issue, one has also take into 
account weight limits, forbidden ways, obligated passages for some categories of 
vehicles, but also many psychological factors such the driver predilection for 
large, high-capacity freeways, the attraction exercised by some areas of the city 
for tourists, seasonal loading of the roads etc. 
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On the other hand, beyond accuracy, a useful traffic model should also be 
evolutionary. Nowadays, the city skyline is continuously changing. Even if large 
freeways do not change for long periods of time, different restrictions based on 
weight and category may be made for vehicles (by example, reserving particular 
lanes for bus traffic only). Also, infrastructure maintenance interventions, natural 
disaster, all require reconfiguring the traffic model, the time being a critical issue. 


A flexible and appropriate model is useful not only for understanding the city 
fluxes of vehicles and analyzing the apparition of bottlenecks and of heavy-traffic 
areas. Ideally, it should also allow the legal authorities to properly react in order to 
minimize the negative (for traffic fluency) effects of unexpected incidents, such as 
collisions, ambulance, police or firemen interventions, or simply to help flushing 
the formed. bottlenecks. This issue is even more important as the congestions 
caused by the above-mentioned usually have side effects, sometimes reverberating 
into the whole traffic. Although the basic causes are usually quite simple (by 
example, one driver may just attempt to use a different route or to turn around), 
the outcome may be a domino-like effect. These interdependencies and resulting 
difficulties mean that the separate modeling of different areas of the city is no 
longer appropriated. 


The unpredictability of city traffic, directly caused by the unpredictability of 
drivers’ intentions and behavior, called for the use of various techniques, ranging 
from statistics to hydrodynamics and mean field approaches [1-8] or chaotic 
models [9] [10]. The latter, by example, is a promising tool to rapidly assert, even 
if with an underlying probability, the effects of events such as car collisions, 
blockages etc. at the level of the entire city. 


2. Peculiarities of city traffic modeling 


Regarding city traffic modeling, an essential requirement is to be able to produce 
a dynamically correct representation of the transportation system. A simple 
average of streets loading will not provide a valuable measure of peak-hours 
traffic and will fail to predict congestions caused by such general hurries as going 
to work and coming back home. Thus, dynamical correctness, along with wide 
temporal and spatial scales is required for the traffic model. Such are the 
microscopic approaches, based on the individual description of the cars, seen as 
the smallest particles in traffic, unlike the method of aggregated traffic flow. 


The microscopic approach has a number of advantages, such as allowing 
including individual route choice behavior for each car. However, even if 
individual car models are rather simple, they quickly sum up, so that modeling of 
large transportation systems [11] may prove computationally very intensive. 


An intermediate way is to provide detailed, i.e. microscopic modeling for critical 
traffic points only. Such are the transitions between the roads of the city. The 
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dynamics of these transitions is important, particularly since it is at the level of 
these transitions that the bottlenecks and traffic congestions tend to occur. 


In urban systems, probably the most important bottlenecks are the traffic lights. 
Usually, these are studied via queuing-type models [12-14]. In order to account 
for the road capacity, a limited service rate is introduced for each link. When the 
time necessary for a car to arrive at the end of the link is reached, the car is added 
to the queue in front of the traffic light. 


A different approach, presented in [15] and reviewed here, makes use of simple 
single-lane micro-simulations in order to capture the dynamics of the traffic on the 
link itself. 


A simulation model of urban traffic is presented below, using a combination of 
stochastic cellular automata and probabilistic transitions between streets. The 
urban road network is represented via an associated graph, where a link means a 
directed street segment (a two-direction road will be divided in two such oriented 
links) and a node models an intersection. Each link of the graph will be bounded 
by an input node and an output node. Vehicles are moved on the links according 
to an underlying cellular automaton model described below, and are transferred 
between two paired links (i.e. having the same input and output nodes) 
probabilistically, based in the link’s capacity. 


3. Models for links and transitions 


Roads (and corresponding links) have different characteristics, modeling their 
physical parameters: length, speed-limit, number of lanes, maximum capacity etc. 
In the approach presented here, from the length of the road and the speed limit 
follows the time necessary for the car to arrive at the end of the road. But, the 
road’s length will also determine the number of the sites (cells of the cellular 
automaton) allowable for the cars on the given road. 


Similarly to [16], the length of one site is assumed to be 7.5 m. Each such site 
(cell) can be either empty or occupied by a vehicle. The latter has a state defined 


by its quantified velocity v € {0...Vina }- A typical choice [15] is to use 5 


quantification levels, so thatv,,,. =5. 


Recall that each link is, in fact, a one-lane segment. Based on this assumption, it is 
assumed that the vehicles are moving using the Nagel-Schreckenberg model [16]. 
The number of unoccupied sites in front of a vehicle is denoted by gap. Normally, 
a large gap will encourage the driver to accelerate. However, there is also the 
possibility for this expected behavior to no occur. This case, i.e. having a driver 
slower than he could, is modeled by the probability Dyvise. 


The model (and its simulation) is iterative. Let us denote by rand a normalized 
(i.e. between 0 andl) random number, generated at each iteration. One such 
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iteration of the traffic model consists of the following three sequential steps, 
applied in parallel to all cars: 

1. all free (having a strictly positive gap) vehicles are accelerated at 
the next level of speed, provided that the maximum allowed speed 
is not already reached: if(v<Vmax) then v=v+1; 

2. the effect of the other cars is taken into account by slowing all the 
cars which are too close of theirs forerunners: if(v>gap) then 
v=gap; 

3. driver decelerating behavior is stochastically modeled: if (v>0) and 
(rand<ppoise) then v=v-1. 

For each link, the transition probability is a function of the capacity of the link. 
This is useful in order to model multi-lane roads (where a road is seen as a set of 
paired links). Obviously, the one-lane model is faster and easier to implement and 
to simulate. 

Differentiating the roads in a city, and the corresponding graph links, is achieved 
by associating probabilities. They help to differentiate the existing links within a 
city, which vary from high capacity segments such as freeways to low capacity 
segments such as arterials. If only one-lane links are considered, the probabilistic 
transition is introduced to control the output flow of a link. A high-capacity link 
will produce a high output flow, while a low capacity link will produce a low 
output flow. 


4. Simulation of a simple node 


In this part, two consecutive links separated by a probabilistic transition Pyrans are 
considered (see Fig. 1). Each n iterations, a vehicle is introduced in the first site of 
the link 1, unless it is already occupied. The iterative model is ran and, when 
vehicles reach the separation between the two links, they are transferred in the 
first site of the second link (only if that site is free), with the probability Dyans. 


Link | Link 2 


Entry point | | | Exit point 


p_trans 
Fig. 1. Two links separated by a transition probability 


In the particular case when the Pyrans equals 0.5, an example of the simulated flow 
(in vehicles per iteration) at the output of the second link is shown in Fig. 2. In the 
simulated scenario, a new car is introduced every 3 iterations in the first cell of the 
first link and a number of 5 speed levels are assumed. If the transition of the 
vehicle towards the second link is not possible at the current iteration (which 
happens with the probability 1- Pirans), then the situation is modeled as the traffic 
light is red. For the vehicle advancing towards separation, the situation is identical 
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as if the first site of the second link is occupied. As such, the iteration 2 of the 
model (i.e. truncating the car’s speed at gap level) will be applied to the vehicles 
having red light. One should be aware that the traffic light color is subject to 
change at every iteration, so the same vehicle may successively accelerate and 
decelerate as approaching the separation. 


output flow (p_noise=0.5 p_trans=0.5) 
0.25 : 
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Fig. 2. Flow at the output of the second link 


The fact that the transition of the vehicle reaching the intersection is not 
guaranteed for all iterations means blockages could occur. This is clearly seen 
with the help of space-time diagrams. These are graphics representing the 
temporal evolution of the underlying cellular automaton. Black points designate 
the occupied cells while white areas mean unoccupied cells. A state-space 
diagram has the spatial information represented by the upwards-oriented vertical 
axis, while the rightwards-oriented horizontal axis represents the time bins. 


Two cases are illustrated: for Pjrans=0.5 (see Fig. 3) and for pyans=0.9 (see 
Fig. 4). 
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Fig. 3. Space-time diagram (Pj,qns=0.5) 
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Simulations show that the output flow measured at the end of the second link and 
the probability Pyans are not proportional. This is because of the non-linearity of 
the assumed model, respectively of the cellular automaton. 
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Fig. 4. Space-time diagram (pjrans=0.9) 


The illustration of this phenomenon is given in Fig. 5, where the average flow is 
depicted versus the values of Pjrans. The plot may be divided into three areas, 
analyzed below. 
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Fig. 5. Average output flow as a function of Dyans 


The first area corresponds to values of Dirans between 0.8 and 1. For such high 
transition probabilities, the output flow varies linearly with pyans. This is expected 
since in this case vehicles do not stop at the separation except very rarely. The 
space-time diagram in Fig. 4 corresponds to this area of the graphic (Pyrans 1S 0.9). 


Cellular automaton urban traffic model 119 


The second area is the leftmost one, corresponding to values of Dirans between 0 
and 0.4. Because of the low transition probability, the vast majority of cars stop at 
the separation. This generates compact traffic jams, such as the ones shown in Fig. 
6 are. 
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Fig. 6. Space-time diagram (Pjrans=0.2). 


It may be shown [15] that, in this area, the corresponding flow is 
P noise — P trans 


ar W? ans 
output flow as long as Pjrans 18 inferior to 0.4. 


approximately F = . This function (also shown in Fig. 5) fits well the 


Finally, the third area is corresponds to values of Pian; between 0.4 and 0.8. In this 
case, an analytical approximation is not as easy to formulate. It is worth noting 
that the flow may be superior to the analytical approximation of proportionality 
(which is valid for large Pyrans values), as shown in Fig. 5. 


5. Simulation of a classic traffic light 


While the previous results allow one to grasp the underlying concepts, it is 
unusual in practice to rely on a randomly-changing traffic light, possibly at every 
iteration. A more realistic situation is to consider a traffic light which keeps the 
fraction of a green light of the total time of a traffic cycle constant. 


From a stochastic point of view, this is analogous with the already-studied case of 
a random traffic light of transition probability Pian; if the approximation 
Prrans=T gree T green+T rea) is used. However, in this case the transition is no more 
probabilistic, but deterministic. Further, the already discussed two-links model is 
analyzed from this perspective. 
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Classical traffic light 


T green Trea T green 


Fig. 7. Transitions of a classical traffic light 
The obtained space-time diagrams are illustrated in Fig. 8 (for Pyrans=0.5) and in 
Fig. 9 (for pirans=0.9), respectively. 


state-space (p_noise=0.5 p_trans=0.5) 
ZT T a: 


TE 


2505 


cell index 


100/- 


50 100 150 200 25 300 35 400 
iteration 


Fig. 8. Space-time diagram for classical traffic light (Dj-4);=0.5) 
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When Pjrans=0.5, the output vehicle flow is quite syncopate and compact traffic 
jams appear. An easy analytical approximation is unlikely. For high values of the 
transition probability, such as Pyans=0.9, vehicles still have to stop occasionally, 
but the output flow is, as expected, superior to the case Pjrans=0.5. 
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Fig. 9. Space-time diagram for classical traffic light (Dj-an;=0.9) 
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However, because of these stops, the output flow remains inferior to the first- 
studied, random case, even for large Pian; values. The analysis of the output flow 
with respect to the Prrans variable, illustrated in Fig. 10, shows an almost linear 
behavior. As already remarked, in the rightmost region the values are just under 
below the linear approximation, while in the random traffic light case, they are a 
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Fig. 10. Average output flow as a function of Djrans 


Conclusions 


Although it represents only a basic tool, the presented approach is useful in that it 
may be quickly extended to complex city road networks. It provides a realistic 
modeling tool for the main blockage sites in a city, i.e. the traffic lights. Such 
model is a valuable tool when it comes to increase fluency of traffic, while 
minimizing fuel consumption and pollution. The approach allows extracting 
analytical approximations, which may be used to speed up computations. 

The simplicity of the model makes it appropriate for large-scale simulations or 
when different alternative traffic configurations are tested and evaluated. It can be 
extended to more complicated roads networks by considering a more elaborate 
transition procedure. 

Also, the proposed model may be used for evaluating the impact on existing 
structures of the projected modifications (new bridges, roads etc.). 
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